Reconstructing a pure state of a spin s through three Stern-Gerlach measurements 
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Consider a spin s prepared in a pure state. It is shown that, generically, the moduli of the 
(2s + 1) spin components along three directions in space determine the state unambigously. These 
probabilities are accessible experimentally by means of a standard Stern-Gerlach apparatus. To 
reconstruct a pure state is therefore possible on the basis of 3(2s + 1) measured intensities. 
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The reconstruction of a particle density-operator is 
possible in principle through repeated measurements on 
an ensemble of identically prepared systems |l],|^] . Quan- 
tum states of vibrating molecules || , of trapped ions j| , 
as well as the state of atoms in motion j5j have been 
reconstructed successfully in the laboratory. Similarly, 
quantum optical experiments |(| have been performed. 

For a spin of length s, this question arises for states in 
a Hilbert space of finite dimension. There is an explicit 
expression for the density matrix p in terms of the moduli 
of spin components along (4s + 1) appropriate directions 
in space |?]]. This number can be reduced to (2s+l) upon 
adopting a different approach ||. A standard Stern- 
Gerlach apparatus with variable orientation in space pro- 
vides the corresponding probabilities in an experiment. 
Alternatively, a Wigner function defined on the discrete 
phase space associated with a finite-dimensional Hilbert 
space allows one to reconstruct quantum states || . This 
method has been adapted in ]To| in order to determine a 
quantized electromagnetic mode of a cavity. Every pro- 
posed method of state reconstruction is bound to reflect 
on the link between the outcomes of a finite number of 
measurements obtained in an actual experiment and the 
mathematical probabilities which refer to infinite ensem- 
bles (see jllj, for example). 

Suppose now that the spin state to be reconstructed is 
known to be prepared in a pure state which is determined 
by less parameters than a mixed one. How to exploit this 
additional knowledge in the most efficient way? Recon- 
struction of pure states has been turned into a question 
as early as 1933 for a particle by Pauli Jl2| who did not 
provide an answer. One solution of the spin version of 
the problem |rq ] makes use of a Feynman filter. This is 
an advanced version of a Stern-Gerlach apparatus which 
is assumed to reveal the relative phases of the expansion 
coefficients of a pure spin state. Another approach relates 
expectation values of spin multipoles with the parameters 
which define the quantum state jlj} . 

As shown in this letter, the pure state of a spin s is 
determined unambigously if the intensities of the spin 
components are measured along three axes. Compared 
to the (2s + 1) axes required for a mixed state the ex- 
perimental effort to perform state reconstruction is thus 
reduced considerably for large spins. Further, this result 



is satisfactory from a mathematical point of view since it 
generalizes an earlier result: the intensities along two in- 
finitesimally close axes spanning a plane define a unique 
pure state when complemented by the expectation value 
of a spin component "out of plane" [Gl| . Effectively, this 
means to measure (2s + 1) probabilities along a third 
direction. 

The states of a spin of magnitude s live in a Hilbert 
space H s of complex dimension (2s + 1), which carries 
an irreducible representation of the group SU(2). The 
components of the spin operator S = hs with standard 
commutation relations [s^,s a ] = is z ,... generate rota- 
tions about the corresponding axes. The standard basis 
of the space HP is given by the eigenvectors of the z com- 
ponent of the spin, denoted by |s,/i z ), — s < ji z < s. The 
transformation under the anti-unitary time reversal op- 



erator T fixes their phases, T\s,/i z ) = (— l) s 
When expanded in the z basis (fik = f-z), 



x,y,z. 
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a pure state is seen to be determined by (2s + 1) complex 
coefficients tp^ = (s,Pz\ip)- If normalized, rays \ip) de- 
pend on 4s real parameters. Two other bases of the space 
H s are used in Eq. (|l]): the sets {\s,fi x )} and {\s,[i y }} 
with — s < fj, x , (j, y < s, made up from the eigenvectors of 
the spin components s x and s y , respectively. Rotations 
about appropriate axes by an angle 7r/2 map them to the 
z basis: 



-lTTSy/2 



8, fix) 



e™ s */ 2 \s,n y ) 



(2) 



A measurement of the intensities {|(s, p z \ip)\ 2 } does 
not fix a single state since the phases of the coeffi- 
cients ipn x remain undetermined. However: 
a spin state \ip) G H s is determined unambiguously if 
3(2s + 1) probabilities 



p(Mfc) = l^wl 2 > k = x,y,z. 



(3) 



are measured with a Stern-Gerlach apparatus along three 
axes not in a plane. For some exceptional states of 
measure zero in Hilbert space H s , the probabilities p(fJ-k) 
might be compatible with a finite number of states. 
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For simplicity, the proof is carried out for orthogonal 
axes, the generalization being straightforward. Measur- 
ing with respect to two axes provides 2 (2s + 1) intensi- 
ties which are usually compatible with a huge number 
of isolated states, in agreement with the result of [ p^[ : 
the parameters fullfil nonlinear relations which may have 
multiple solutions. Enumerating the ensemble of possi- 
ble "partner" states is complicated, so a distinctive third 
measurement is included from the very beginning. 

It is useful to rephrase the statement at stake differ- 
ently According to (|J) a state \ip) gives rise to the same 

intensities as does |V>) if its coefficients V/j fc = (s i A*fe|V') 
differ from by phase factors only. Using (|l|) one 
writes thus 



E 



Vwe <Xfe(Wt) |s,Mfc> = exp[z Xfc ( Sfc )]|V) , (4) 



with three polynomials Xfe(A*) °f order 2s in /x at most. 
From now on, the index k is understood to take the values 
x, y, and z throughout. The coefficients in (^) thus define 
three states \ipk) = W% \ip), where Wjl = exp [ixk( s k)] is 
a unitary operator diagonal in the k basis. Consequently, 
a state \ip) compatible with (||) exists if and only if there 
are nontrivial unitary operators Wt such that 



w°\i>) = w°\i>) = w°\i>) = 



(5) 



It will turn out that this relation is satisfied only if the 
operators W£ are multiples of the identity, implying that 
\tp) and \ip) represent the same ray in Hilbert space. 

Before turning to the proof, the intensities p((J-k) in 
(||) are represented in a more compact way. Define three 
functions mk(a) of a complex variable aeCby 



m k (a) = (iJj\U^a)\4>) 



(6) 



where the operator Ujl(a) — exp(iasfc) rotates a state 
\ip) about the A; axis if a € R. Eq. (§) is inverted easily 
using the orthogonality of the functions exp[— i/x^a] on 
the interval < a < 2tt. 

The proof showing that the data (||) are sufficient for 
state reconstruction is divided into five steps, (i) A 2 2s 
dimensional "parent" space TL S is introduced which con- 
tains the Hilbert space H s of the spin s as a subspace. (n) 
To each state V) G H s an equivalence class of product 
states {\ty) G H s } is associated, (m) A natural definition 
of generic states emerges for product states in TL S and, 
a fortiori, in H s . (iv) An appropriate set of expecta- 
tion values of the parent states fixes them uniquely, 
(i?) Finally, it is shown that all states \ip) satisfying (|^) 
have parents in the same equivalence class as the original 
IV')- Consequently, the (generic) state \ip) is the only one 
giving rise to the intensities (0). 



(i) The 2 2s dimensional "parent" space H s of H s is 
obtained from tensoring 2s copies of the Hilbert space 
C 2 of a spin 1/2: 
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Cl 



(7) 



r=l 



A basis of C 2 is given by the eigenstates \<r) = \s = 
1/2, /i3 = a/2), a = ±1, of the third component of the 
spin 1/2: er 3 |er) = <j\<j). This choice induces a basis of 
Ti s formed by all product states 



2 s 

IK»=(g) Wr). 
r=l 



(8) 



The parent space H s decomposes into a subspace 7Yf ym 
and its complement, 



(9) 



where 7if ym is spanned by the (2s + 1) states obtained 
from completely symmetrizing those in (^): 

|s,M3) = S 2s \Wr}) 

S{ai + --- + °~2s -2^ 3 )|{CT r }), (10) 

{a r } 



Nf. 



where — s < /i 3 < s, using a symmetrizer of 2s objects, 
S2 S , and the normalization factor iV* = ((s — /i 3 )!(s + 
/ u 3 )!/(2s)!) 1 / 2 . The space 7Yf ym is important here because 
it carries a (2s + l) dimensional irreducible representation 
of the group of rotations, SU (2) , obtained upon reducing 
the product representation 16 



U\{* r }) 
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(11) 



=1 <j'=±i 



where u r is the r-th copy of a rotation u G SU (2) of the 
fundamental representation acting on C 2 , and U is an 
operator defined on H s . Since Hilbert spaces of the same 
dimension are isomorphic, W| ym and H s will be identified 
from now on. 

(n) There is a one-to-one relation between states |V) G 
7if ym and equivalence classes of product states |^) G H s : 



Is / \ 

\#) = \{* r }) = ®[Y l KMr) 
r=l \ oy / 



(12) 



The equivalence relation ~ is defined as follows: the 
projection of a state \ty) in ( |l2| ) onto a basis state 
s,/i 3 ) G Htym must equal the corresponding expansion 
coefficient of V) in the z basis, i.e. , 



(s,/i 3 |*) = AV(s,m z |V>) 



-s < Hz = Hz < s . 



(13) 



and the factor iV^ > may depend on the state V) under 
consideration but not on the index fj, z . Thus, ~ |^') 
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means that for a fixed |?/>), the Eqs. ( J13| ) hold for both 
product states, ~ |\E''). The association of spin states 
with product or "parent" states \*f?) is essential for 
the following. 

In order to determine the class of states satisfying Eq. 
( |l3| ) for a prescribed vector (with definite phase), 
multiply by the factor 1/iV*, by powers (— z) M+s and sum 
all terms. The right-hand-side then defines an analytic 
function 



f R (z) = J2 



(-*) 



ft+S 



2s 



fl — — S 



OC JJ(z r - z) . 



(14) 



specified by the location of its 2s zeroes z r in the com- 
plex plane. The left-hand-side yields a second analytic 
function of z, 



fl=-s {ay} 
2s 



vl/ 2s 



(15) 



one sees that the ensemble {z r } = (z%, Z2 S ) (no or- 
der implied) defines the daughter completely while a 
maximum of (2s)! different parent states |\&) is associated 
with a given set {z r }. 

(in) Suppose that three ensembles of 2s real num- 
bers each, {x r }, {y r }, and {\z r \} = (\zi\, \zzs\) with 
z r = x r + iyr are given in disorder. If one is able to 
construct the disordered ensemble of 2s complex num- 
bers {z r = x r + iy r } upon using the 2s conditions 
\z r \ 2 = x 2 + y 2 , the equivalence class with represen- 
tant \^f) is called generic. In other words, it must be 
possible to combine unambiguously real and imaginary 
parts into complex numbers z r . In this spirit, a daughter 
\ip) G W| ym will be called generic if it has generic par- 
ents The procedure does not work if equalities such 
as x r = ±y r ',r ^ r' exist; hence exceptional states have 
measure zero. 

(iv) It is shown now that the expectation values of rota- 
tions Uk(ct) about the axes x, y, and z, fix generic product 
states |^) = |{vE ,r }) up to a permutation of the factors 
\^ r ) and an overall phase factor. A generic \^f) 6 H s 
leads to three expectation values 



The (2s + 1) equations ( |l5| ) are satisfied if fh(z) and 
/r(z) coincide. Being two polynomials of degree 2s, this 
requires them to have identical zeroes, 



,2s: 



(16) 



in addition, /z,(0) = /_r(0) must hold. Due to the nor- 
malization (^ r \W) = I 2 + |\F I 2 = 1, one can write 



^l + \z r \ 2 \z r 



1 



sr G [0, 2tt) . (17) 



M fc (a) = (*|W fc (a)|*> = l[(* r \u r k (aW 



(20) 



r=l 



where Uk(a) = lcos(a/2) + er^ sin(a/2) represents a ro- 
tation about axis k in C 2 . Using the parametrization of 
Eq. (p^), the functions Mk(a) defined in ( p0| ) read ex- 
plicitly 



M x (a) = Yl 



cos(a/2) + 2ix r sin(a/2) 

l+\z r \ 2 



(21a) 



Thus, there are 2s undetermined phase factors e lKr with a 
product equal to 1 (remember that \tf>) denotes a vector). 
However, the overall ambiguity is even larger: when com- 
paring the zeroes of the functions Jl(z) and fn(z), there 
is no rule which would indicate what order to choose 
when writing down the product state {'I' 1 *}). In other 
words, the equivalence class of states defined by (|l^) con- 
sists of all states with coefficients ([l7|) distributed in any 
order over the 2s spinors in (|l2[). All these states are 



parents of the same since they satisfy Eq. (13) 
A given product state |\&) with components 



<K}I*> = * 



(18) 



2.s 



M y (a) = H 



cos(a/2) + 2y r sin(a/2) 



r=l 



l + |z,f 



(21b) 



£j cos(a/2) + i(l - \z r \ 2 ) sin(a/2) 

M z {a) = [[ rro ' ( 21c ) 



1 + Zr 



where again z r — x r + iy r . Denote by = |{\P r }) 
another product state with expectations Mk(a): 



2.s 



M k (a) = (*|W fc (a)|*> = H$ r \K(<xW) ■ (22) 



has a unique "daughter" \ip) to be read off directly. Upon 
parametrizing each factor \^ r ) by a complex number z r , 



*+ 1 = 1 { 1 



(19) 



Upon describing the state \ i $>) by the sequence {z r }, the 
three functions Mfc(a) are given by Eqs. |l|) after replac- 
ing each z r by z r . It is shown now that the conditions 

(*|W fc (a)|*> = (*|Wfc(a)|*>, k = x,y,z. (23) 
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necessitate |^) ~ |\&). Being analytic in the complex 
a plane, the functions M^(a) and Mk(a) are equal if 
they have same zeroes. The equality M z (a) — M z (a) 
requires \z r \ = \z r \. The condition M x (a) = M x (a) 
in^turn implies x r — x r ; finally, y r = y r follows from 
My (a) = M y (a). However, this procedure determines 
the ensembles {x r }, {y r }, and { 1 2r r | } without any order of 
its members. Nevertheless, one can reconstruct the en- 
semble {z r } (no order implied) according to (in) if |\&) 
is generic providing thus a unique equivalence class. For 
exceptional states, the 2s complex numbers cannot be 
reconstructed unambigously since they might allow for 
parents contained in different equivalence classes. 

(v) The results (?) to (iv) imply that the probabilities 
p(pt,k) for three directions x, y, and z as given in Eq. (j^) 
determine a generic state \ip) unambigously. According 
to Eq. (||), a state \ip) gives rise to the same probabilities 
as does \ip) if one has 

|<M = = \1>.) = \$) ■ (24) 
For parent states \&k) of \ipk) this relation says that 



(25) 



This implies that the mean values {^! k\lA x (a)\%! k) of 
the operator U x (a) = ® r exp[iaa x /2\ are equal for k = 
x,y,z: as products they are invariant under a permu- 
tation of their factors. This also holds for expectation 
values of the operators U y (a) and U z (a). Write the 
parent states Y&k) in the form Wfc|vE') with operators 
yVk({&k,r}) = <8V exp[iafe !r cr J i £ /2] defined on the parent 
space 7i s such that they have W% as component acting 
in 7i| ym . Contrary to the rotations lAk(a) which depend 
linearly on the generators s^, the operators Wi are non- 
linear functions x(s^) of them, Eq. (Q). Therefore, the 
operators Wk({ctk,r}) depend on a set of 2s different an- 
gles {ctk, r }- Using (|25|) one concludes 



(*|«*|*> = (* fc |Z4|*fc) 

= (v\wtu k w k \*) = (*|W t |*) 



(26) 



The third equality follows because Wfc and Uk do com- 
mute, both being functions of 8k only. Eq. (pgj ) comes 
down to saying that the functions Mfc(a) and M k (a) co- 
incide for all k and a. One concludes thus with (w) that 
the state \^), a parent of is necessarily a member 
of the same equivalence class as the parent \^f) of 
In other words, the application of the operators Wfc on 
a parent |\&) does not map it into another equivalence 
class. In the generic case, there is thus no state different 
from \ip) with the same data (||) what was to be shown. 

The reasoning (i) to (v) remains valid if one measures 
the intensities along directions characterized by unit vec- 
tors 11^,11,,, and instead of three orthogonal axes. 



These vectors must be linearly independent, that is, they 
have to span a volume in space: • x =/= 0. 

As a matter of fact, it is not excluded that the set 
of data (||) be also sufficient to determine exceptional 
states unambiguously. Suppose that the numbers {z r } 
are associated with a parent state l^) and {z' r } with 
another one, where both sets of complex numbers 

are obtained from the ensembles {x r } and {y r } through 



\Zr\< 



yf.. This does not necessarily imply the ex- 



istence of an independent ^ \ip) since it is the basic 
conditions |V>' | = \ip^ k | which must be satisfied. Explicit 
calculations for low values of spin s show that this hap- 
pens only if 7// fe = tp* h , resulting in {ip\(s y ) 2n+1 \i}}) = 
for all integers n. In any case one expects every non- 
genericity to vanish if the spatial directions involved are 
slightly modified. 

To sum up, state reconstruction is possible if based on 
the 3 (2s + 1) moduli of the spin components with respect 
to three directions in space not all in the same plane. 
Compared to a constructive method using (2s + l) 2 real 
numbers, the non-constructive method presented here re- 
quires that considerably less parameters be determined 
experimentally, namely 3(2s + 1). 
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